0D£  m COPY  , . ADA  0383  0 6 


DYNAMIC  FRACTURE  AND  FRAG  M E NT  AT  I ON 
OF  CYLINDERS 


by 

C.  ROBERT  CROWE 
WILLIAM  H.  HOLT 


j 


NAVAL  SURFACE  WEAPONS  CENTER 
DAHLGREN  LABORATORY 
Dahlgren,  Virginia 
22448 


D.  M.A*new,  Jr.,  Capt.,  USN  j H(  MiU5|  j, 

OIC  and  Aaslalant  Commander  Associate  Technical  Director 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whir,  Polo  Entered) 

REPORT  DOCUMENTATION  PAGE  befoI^compleIinTform 

1.  R€PO«T  NUMBER  (^GOVT  ACCESSION  NO-  3*  RECIPIENT’S^ AT ALOG~NUMBER 


1.  R€PO«T  NUMBER 


TR-3449 

4 TITLE  (and  Subtitle) 


5.„  TyPE  OF  REPORT  ft  PERIOO  COVERED 


1 I 


DYNAMIC  FRACTURE  AND  FRAGMENTATION  OF  CYLINDERS , \ Final  -r  < y't  * 


It.  nuTwoRr*; 


6.  CONTRACT  OR  GRANT  NUMBERfaJ 


B’PCRFiORMIRG  ORGANIZATION  NAME  AND  ADDRESS 

Naval  Surface  Weapons  Center^ 
Dahlgren  Laboratory 
Dahlgren,  Virginia  22448 

II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  A WORK  UNIT  NUMBERS 


\IZ.  REPORLUBIE — 

f November  M76  j 


HUMBER  OF  PAGES 

65 

U.  MONITORING  AGENCY  NAME  A AOORESSflf  dUtoront  Iraa i Conlrolllnl  officbj  IS.  SECURITY  CLASS,  (ol  Ihle  report) 


/>  NS  vc/s  A - T H - Tvy? 

'i — ■ ■ ■■  ~ 

16.  DISTRIBUTION  STATEMENT  (at  Ihle  Raporl)  — 


UNCLASSIFIED 

13a.  DECLASSIFICATION/DOWNGRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (of  thle  /import)  t ] 

( /> J (Si//--* 

Approved  for  public  rejwa^e;  distribution  unlimited.  / 


Approved  tor  public  reJ^ase;  distn 

m.i-ijL S'Sfe-. 


* 553  5 p -■-{ 


I 17.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  In  Block  20,  II  dlfforont  from  Roport) 


SUPPLEMENTARY  NOTES 


19*  KEY  WOROS  (Continue  on  tovoroo  oldo  If  neceeemry  and  Identity  by  block  nuai bar) 

Fracture  Solid  Mechanics 

Fragmentation 

Shock  Waves 

Shear  Banding 


80.  APS!  RACT  (Continue  on  reeeree  old e II  nniil»T  ond  Idonllty  br  block  number) 

An  experimental  and  theoretical  investigation  of  the  mechanisms  of  fragmenta- 
tion has  been  conducted  for  explosive-filled  cylinders  of  Armco  iron  and  HF-1 
steel.  Dynamic  fracture  nucleation  and  growth  parameters  for  HF-1  steel  were 
determined  from  gas  gun  experiments  for  use  in  computer  modeling  of  the  frac- 
ture and  fragmentation  processes.  Calculated  fragment  mass  distributions  are 
compared  with  experimental  results. 

\ 


DO  1473  COITION  OF  1 NOV  69  IS  OBSOLETE 

S/N  0102*  LF-  014-  6601 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  fW»«n  Data  Entered) 


/ • 


FOREWORD 


A joint  NSWC/DL  and  Stanford  Research  Institute  (SRI)  experimental 
and  computational  effort  has  been  carried  out  to  develop  a quantitative 
description  of  the  fracture  and  fragmentation  of  explosive-filled  cylin- 
ders. The  approach  is  based  on  previously  developed  SRI  models  for 
fracture  processes.  This  report  presents  the  NSWC/DL  part  of  the  work 
and  briefly  describes  the  SRI  part.  The  work  was  supported  by  NAVSEA 
Task  No.  SF  32-353-302  and  by  NSWC/DL  Independent  Research  Funds. 

This  report  has  been  reviewed  and  approved  by  C.  A.  Cooper,  Head, 
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EXECUTIVE  SUMMARY 


A joint  NSWC/DL  and  Stanford  Research  Institute  (SRI)  experimental 
and  analytical  effort  has  been  carried  out  to  predict  the  fragmentation 
performance  of  explosive-filled  metal  cylinders.  The  approach  is  based 
on  computational  models  that  have  been  developed  by  SRI  to  describe  the 
radial  fracture  and  shear  deformation  processes  that  occur  in  cylinder 
fragmentation.  The  approach  utilizes  very  high  strain  rate  material 
parameters  that  are  determined  from  gas  gun  impact  experiments  on  small 
specimens.  These  parameters  are  used  in  computational  models  for  the 
nucleation  and  growth  of  cracks  in  the  material.  Thermo-mechanical  data 
are  used  in  modeling  of  the  shear  deformation.  The  procedure  has  been 
applied  to  explosive-filled  cylinders  of  Armco  iron  and  HF-1  steel. 
Computational  results  for  the  fragment  mass  distributions  show  favorable 
agreement  with  experimental  distributions. 
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The  design  of  a fragmentation  warhead  is  influenced  by  many  factors, 
including  cost,  logistics,  safety,  and  effectiveness  in  damaging  a tar- 
get. The  design  objective  is  to  maximize  target  damage,  consistent  with 
the  other  factors.  A basic  problem  is  to  design  the  warhead  to  have 
fragments  with  the  mass  distribution,  velocities,  and  angular  distri- 
butions needed  to  defeat  the  target  in  an  anticipated  type  of  encounter. 
Experience  derived  from  years  of  extensive  testing  of  prototype  warheads 
has  shown  that  the  fragment  directions  are  controlled  mainly  by  warhead 
geometry;  however,  the  fragment  mass  distribution  depends  on  the  metal 
case  thickness,  explosive-mass  to  metal-mass  ratio,  and  the  metallurgi- 
cal properties  of  the  case  material.  While  these  general  facts  are 
known,  adequate  theoretical  models  for  predicting  the  fragmentation 
behavior  of  a given  warhead  configuration  have  not  evolved.  A predictive 
capability  for  fragmentation  would  permit  the  computational  optimization 
of  warhead  design  for  given  target  and  launch  conditions,  and  would 
minimize  the  costly,  time-consuming  testing  of  prototypes. 

A new  quantitative  approach  to  the  description  of  fracture  and 
fragmentation  processes  has  led  to  a predictive  capability  for  armor 
penetration  problems;'*’  the  work  described  herein  is  part  of  a joint 
effort  by  NSWC/DL  and  the  Stanford  Research  Institute  (SRI)  to  apply 
this  approach  to  predict  the  fragment  mass  distribution  (Mott  plot)  for 
an  explosive-filled  metal  cylinder.  The  new  approach  involves  experi- 
mental determination  of  the  very  high  strain  rate  (~10'’s  *)  fracture 
response  of  a material  and  computational  modelina  of  the  fracture  and 
fragmentation  processes.  Fracture  data  for  a material  are  obtained 
from  gas  gun  impact  and  soft  recovery  experiments  on  small  specimens. 

These  data  are  used  to  obtain  crack  nucleation  and  growth  (NAG)  parameters 
that  quantify  the  processes  leadinq  to  fragmentation. 

When  the  explosive  in  a cylindrical  warhead  is  detonated,  the  initial 
stresses  in  the  cylinder  wall  are  caused  by  explosive  shock  waves;  these 
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(stresses  cause  rapid  outward  acceleration  of  the  wall.  After  the 

initial  shock  waves  have  subsided,  the  stresses  in  the  cylinder  wall  are 
mainly  those  caused  by  the  outward  displacement.  As  cracks  are  formed 
and  grow,  the  fracture  surfaces  eventually  intersect  to  form  surfaces 
of  individual  fragments.  Previous  investigations  indicate  that  there 
are  two  types  of  macroscopic  fracture  contributing  to  the  fragmentation 
I of  cylinders:  localized  shearing  and  radial  cracking  (Figure  1) . 

RADIAL 


Figure  1.  Two  types  of  fracture  occuring  in  cylinder  fragmentation: 

(a)  pure  shear  fracture  and  (b)  bimodal  fracture. 

As  part  of  the  joint  experimental  and  analytical  effort,  a computa- 
tional model  for  predicting  the  fragmentation  of  a cylinder  was  developed 
by  SRI.  The  brittle  fracture  and  fragmentation  (NAG/FRAG)  model3  pre- 
viously developed  by  SRI  was  modified  to  apply  to  the  radial  cracking. 

A shear  nucleation  and  growth  (SNAG)  model  was  developed  for  the  shear 
fractures.  Details  of  the  SNAG  model  development  are  given  in  Reference 
4.  The  NSWC/DL  work  in  support  of  this  effort  included  gas  gun  experiments 
to  determine  fracture  nucleation  and  growth  parameters  for  HF-1  steel  and 
explosive  fragmentation  experiments  to  obtain  fragment  mass  distributions 
ror  cylinders  of  Armco  iron  and  HF-1  steel. 
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A brief  historical  background  of  cylinder  fragmentation  is  presented 
in  Chapter  IX.  Chapter  III  contains  the  conventional  mechanical  properties 
for  Armco  iron  and  HF-1  steel.  Determination  of  the  HF-1  steel  NAG  para- 
meters is  discussed  in  Chapter  IV.  The  fragmentation  experiments  are 
described  in  Chapter  V.  Shear  band  models  and  a comparison  of  the  experi- 
mental and  SRl-calculated  fragment  mass  distributions  are  given  in 
Chapter  VI.  Chapter  VII  contains  a summary  and  recommendations  for 
future  work. 

II.  HISTORICAL  BACKGROUND 

During  World  War  II,  considerable  research  was  directed  towards 
solving  the  fragmentation  problem.  Most  notable  are  the  works  of 
Taylor^'  6 and  Mott.^  ^ Taylor  devised  models  to  describe  the  fragment- 
ation of  cylindrical  bombs  and  analyzed  the  stress  state  associated  with 
the  cylinder  radius  at  which  fractures  penetrated  the  wall.  He  presented 
data  showing  that  detonation  gases  did  not  escape  through  the  wall  until 
a cylinder  had  expanded  to  nearly  double  its  initial  diameter.  He  showed 
that  the  initial  explosive- induced  circumferential  stresses  near  the  inner 
surface  of  the  cylinder  are  compressive.  Thus,  radial  fractures  that  are 
started  at  the  outer  surface  by  tensile  stresses  could  be  expected  to 
extend  into  the  material  only  to  a depth  at  which  the  stresses  become 
compressive.  Shear  fractures  could  occur  in  the  compressive  region,  but 
these  fractures  would  not  open  to  allow  gases  to  escape  until  the  cir- 
cumferential compression  becomes  zero.  The  predictions  of  this  theory 
are  in  general  agreement  with  experimental  observations,  although  there 
are  some  discrepencies  with  regard  to  details. 

Mott  also  calculated  the  stress  distribution  in  a rapidly  expand- 
ing cylinder  and  used  a Griffith'*'^  crack  model  in  an  attempt  to  predict 
the  propagation  of  cracks  in  the  cylinder  wall.  He  also  developed  a 
theory  for  the  distribution  of  radial  fractures  in  tubular  bombs  that 
were  constructed  of  stacked  rings  around  a thin  cylindrical  shell  and 
filled  with  explosive.  The  theory  incorporates  the  effects  of  local 
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stress  reduction  due  to  relief  waves  that  propagate  from  the  fracture 

surfaces  as  soon  as  a crack  forms.  The  cracks  were  assumed  to  start  at 

inherent  material  flaws  that  had  a Gaussian  size  distribution.  While 
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this  theory  is  generally  assumed  to  be  inadequate,  Mott  and  Linfoot 
found  experimentally  that  the  fragment  mass  distribution  for  most  pro- 
jectiles and  cylindrical  warheads  can  be  approximated  by  an  empirical 
equation  of  the  form 

1/2 

N(>m)  = Nq  exp  [ - (m/y)  ] (1) 

where  N(>m)  is  the  number  of  fragments  with  mass  greater  than  m,  Nq  is 
the  total  number  of  fragments  from  the  warhead,  m is  the  fragment  mass 
group,  and  y is  a fragmentation  size  parameter.  The  exponent  1/2  in 
Equation  (1)  applies  for  fragmentation  processes  where  the  wall  thickness 
controls  one  of  the  fragment  dimensions.  When  the  wall  is  thicker  than 
the  average  fragment  thickness,  the  exponent  1/3  is  used.  Figure  2 is 
a plot  of  Equation  (1) . The  representation  of  fragmentation  data  in 
this  form  has  been  in  use  since  Mott’s  original  work. 


(FRAGMENT  MASS)'/4 


Figure  2.  Graphical  representation  of  the  fragment  mass 
: distribution  from  an  exploded  cylinder 

■ (Mott  plot) . 
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There  have  been  many  other  attempts  to  relate  the  dynamic  processes 
in  fragmentation  to  the  warhead  geometry,  warhead  case  material  properties, 
and  the  properties  of  the  explosive. 

Rinehart  and  Pearson**-  have  proposed  a model  of  fragmentation  in 
which  fracturing  occurs  when  the  radial  particle  velocity  in  the  cylinder 
material  exceeds  a critical  straining  velocity.  The  model  predicts  the 
number  of  fragments  to  be  proportional  to  the  ratio  of  the  radial  particle 
velocity  and  the  critical  straining  velocity. 

A more  thorough  treatment  of  the  fragmenting  bomb  problem  has  been 

12 

carried  out  by  Hoggatt  and  Recht  using  basically  the  same  approach  as 
used  by  Taylor.  They  assumed  that  the  high  pressure  from  the  explosive 
causes  small  radial  cracks  to  appear  in  the  tensile  circumferential  stress 
region  near  the  outer  surface  of  the  cylinder.  Within  this  stress  region, 
cylinder  expansion  occurs  by  extrusion  which  actuates  shear  planes  in  the 
i'.ier  compressive  circumferential  stress  region.  This  model  treats  the 
stress  distribution  in  detail  and  includes  the  effects  of  unstable  thermo- 
plastic shear.  The  unstable  shear  occurs  when  the  material  in  the  shear 
zone  is  thermally  softened  by  the  heat  generated  in  plastic  deformation; 
the  local  flow  stresses  decrease  with  increasing  strain. 

Banks'* ^ has  observed  that  the  number  of  cracks  in  a cylinder  in- 
creases linearly  with  the  radial  velocity  of  the  cylinder  wall  at  the 
time  of  break-up,  as  was  predicted  by  Mott.  In  this  work,  no  considera- 
tion was  given  to  the  absolute  number  of  cracks  and  the  velocities  con- 
sidered were  not  representative  of  those  normally  observed  in  fragmenta- 
tion weapons. 

High-speed  photography  of  fragmenting  cylinder  by  Rinehart  and 

Pearson,*^  Kronman,*-*  and  Clark***  has  shown  that  cracks  start  at  an 

expansion  ratio  of  approximately  1.2  and  detonation  gases  appear  from 

the  cracks  at  an  expansion  ratio  of  approximately  1.8.  These  observations 

are  consistent  with  the  predictions  of  Taylor  and  those  of  Hoggatt  and 
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Recht.  In  a series  of  experiments  by  Clark  and  Juriaco,  it  was  shown 
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that  the  number  of  cracks  in  an  expanding  cylinder  can  be  correlated 
with  Mott's  predictions  to  an  accuracy  of  +20  percent  through  the  use  of 
room  temperature  static  mechanical  properties  and  one  adjustable  para- 
meter. 

18 

Meinert  has  attempted  to  relate  fragmentation  behavior  to  the 

static  fracture  energy,  defined  as  the  area  under  a conventional  stress- 
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strain  curve.  Bardes  reports  that  a linear  relationship  exists  between 
the  y of  Equation  (1)  and  the  percent  reduction  in  area  for  a conventional 
tensile  test  for  HF-1  steel.  Work  by  Magis^’  ^ and  others22  2^  has 
been  directed  toward  correlating  fragment  size  distributions  with  com- 
position, heat  treatment,  and  mechanical  properties.  Magis  considered 
eight  parameters:  charge-mass  to  metal-mass  ratio  (c/m),  outside  dia- 

meter, tensile  strength,  yield  strength,  percent  reduction  of  area,  per- 
cent elongation,  Charpy  impact  strength,  and  hardness  for  10  different 
steels  including  three  heat  treatments  of  HF-1  steel.  He  found  that 
only  three  parameters  (c/m,  outside  diameter,  and  percent  elongation) 
were  statistically  correlated  for  all  the  steels.  The  other  parameters 
had  statistical  correlation  only  in  individual  steels. 

Bardes  investigated  the  effects  of  specific  microconstituents  on 

fragmentation  behavior.  This  work  showed  that  adiabatic  shearing  is 

essentially  independent  of  microstructure,  carbide  networks  at  grain 

boundaries  affect  fragment  size,  large-grained  structures  produce  small 

fragments,  homogenation  produces  large  fragments,  and  decarburization  has 

no  effect.  The  effect  of  the  carbide  network  has  also  been  noted  by 

25 

Clark  and  Juriaco;  and,  work  by  O'Shea  confirms  the  effect  of  grain 
size  on  fragmentation  behavior. 

More  recently,  effort  has  centered  around  predicting  fragmentation 

behavior  through  the  use  of  computer  codes.  The  HEMP  and  FRAG  2 codes 
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have  been  used  in  fragmentation  studies  at  NSWC/DL.  In  HEMP  calculations 
fragmentation  is  introduced  artifically  through  an  additional  operation: 
an  experimentally  determined  average  fragment  mass  for  a Mott  distribution 
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is  applied  to  the  HEMP  results.  This  approach  yields  a distribution  that 
agrees  with  experimental  results  when  fragments  from  all  polar  zones  are 
considered  together,  provided  the  correct  average  mass  has  been  selected 
for  the  computation.  This  procedure  does  not  match  the  distribution  for 
individual  polar  zones.  FRAG  2 utilizes  the  Magis  data  to  calculate  the 
distribution  per  polar  zone  by  calculating  a distribution  for  elemental 
rings  spaced  along  the  length  of  the  projectile,  and  by  utilizing  a 
scaling  equation  to  correct  for  changes  in  ring  size.  The  results 
obtained  with  FRAG  2 are  inadequate  because  of  the  limitations  of  the 
Magis  data  and  the  inability  to  predict  the  experimental  mass  distribu- 
tion per  polar  zone. 


III.  MATERIALS 

The  Armco  iron  used  for  this  work  was  obtained  from  the  Corey 
Steel  Company  as  a 12. 7-cm-diameter  hot-rolled  bar.  The  HF-1  steel  was 
obtained  from  Norris  Industries  as  12. 7-cm-diameter  cylinders  and  hot- 
rolled  bar.  The  manufacturer's  chemical  analysis  of  the  Armco  iron  and 
the  NSWC/DL  chemical  analysis  of  the  HF-1  steel  are  given  in  Table  1. 

Table  1.  Chemical  analysis  of  materials. 


Element 


(Wt  %) 

Armco  Iron 

HF-1  Steel 

c 

0.004 

1.08 

p 

Tr 

0.011 

s 

0.025 

0.0001 

A1 

0.  0002 

0.  006 

Si 

Tr 

0.86 

Mn 

Tr 

1.83 

Mo 

Tr 

0.07 

Ni 

0.0001 

Tr 

Cr 

Tr 

0.15 

0 

0.010 

— 

N 

0.  001 

— 

Cd 

0.003 

— 

Mg 

0.  001 

— 
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The  Armco  iron  was  used  in  the  as-received  condition.  The  structure, 
as  shown  in  Figure  3(a),  is  basically  a single-phase  polycrystalline 
aggregate  with  ASTM  grain  size  number  of  approximately  2 . 5 to  3 . A few 
inclusions  are  also  present;  these  probably  are  oxides.  The  final  structure 
for  the  HF-1  steel  after  heat  treatment  is  shown  in  Figure  3(b).  It 
consists  primarily  of  tempered  martensite  with  some  dispersed  carbide 
precipitates  and  retained  austenite.  Stress-strain  curves  for  Armco 
iron  and  HF-1  steel  are  shown  in  Figures  4 and  5,  respectively.  Relevant 
mechanical  data  for  these  materials  are  listed  in  Table  2. 


Figure  3.  Microstructure  of  materials:  (a)  Armco  iron  and 

(b)  HF-1  steel. 
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Table  2.  Mechanical  properties  of  materials. 


Armco  Iron 
HF-1  Steel 


Ultimate 

0. 2-Percent 

Tensile 

Yield  Stress 

Strength 

(MPa*) 

(MPa) 

Hardness 

180 

300 

Rb  40 

1040 

1320 

R 40 

* 1 megapascal  (MPa)  = 145.0  psi,  1 gigapascal  (GPa) 


10" 


MPa 


IV.  DETERMINATION  OF  FRACTURE  NUCLEATION  AND  GROWTH  PARAMETERS  FOR  HF-1 
STEEL 

A.  INTRODUCTION 

This  chapter  describes  the  method  for  obtaining  the  NAG  parameters 
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for  HF-1  steel.  The  method  developed  by  SRI  for  other  materials  ' has 
been  followed,  with  minor  exceptions.  The  parameters  are  used  in  the  SRI 
brittle  fracture  model  to  predict  radial  crack  formation  in  cylinder 
fragmentation  (Figure  1) . The  steps  to  obtain  the  NAG  parameters  are 
described  briefly  as  follows:  Instrumented  gas  gun  experiments  are 

performed  to  determine  the  shock  wave  equation  of  state  of  the  material. 

A series  of  soft  recovery  gas  gun  experiments  is  performed  to  produce 
differing  amounts  of  internal  fracture  in  the  specimens  under  known 
shock-loading  conditions.  The  recovered  specimens  are  sectioned  on  a 
diameter  and  polished  to  reveal  the  microscopic  spall  fracture.  The 
number  of  cracks,  their  sizes,  and  orientations  with  respect  to  the 
impact  plane  are  digitized.  These  surface  distributions  are  converted 
to  volume  distributions  using  the  SRI  BABS2  computer  code.  The  stress 
state  in  the  material  is  determined  with  the  SRI  PUFF  code.  PUFF  is  a 
finite-difference  wave  propagation  code  that  calculates  the  stress  state 
in  stepwise  time  increments.  A PUFF  computation  is  made  using  the  NAG 
model  and  a first-estimate  set  of  parameters.  This  computation  yields  a 
volume  distribution  of  cracks  which  is  compared  with  an  experimental 
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distribution;  the  values  of  the  NAG  parameters  in  PUFF  are  adjusted  until 
agreement  is  obtained. 


B.  DYNAMIC  FRACTURE  EXPERIMENTS 

28 

Instrumented  gas  gun  experiments  were  performed  to  measure 

the  dynamic  yield  stress  (also  called  the  Hugoniot  elastic  limit)  and 

the  elastic  wave  velocity  in  HF-1  steel.  A Hugoniot  elastic  limit  of 
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2.3  GPa  and  an  elastic  wave  velocity  of  5.94  km/s  were  measured. 

When  performing  the  spall  fracture  experiments,  the  impacted 
specimen  is  soft  recovered  to  minimize  any  unintentional  damage.30  This 
is  accomplished  by  the  post-impact  separation  and  in-flight  capture  of 
the  projectile  end  impactor;  the  specimen  is  soft  recovered  in  a container 
filled  with  rags.  A target  assembly  for  a spall  fracture  experiment  is 
shown  in  Figure  6.  The  gun  bore  diameter  is  40  mm.  The  projectile 
velocity  is  measured  at  impact  with  the  three  velocity  pins.  The  barrel 
is  evacuated  to  minimize  gas  cushion  effects  on  impact.  The  specimen  is 
held  in  a stainless  steel  mounting  ring  by  a low-melting-point  alloy. 


Figure  6.  Schematic  of  muzzle  region  with  specimen  assembly  for 
spall  fracture  experiment. 
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For  spall  fracture  to  occur,  the  experimental  conditions  must 
cause  tensile  stresses  in  the  specimen.  Figure  7 is  a simplified  schematic 
of  the  spall  process.  The  impactor  thickness  is  less  than  the  specimen 
thickness.  In  the  first  step  after  impact,  a plane  shock  wave  propagates 
into  both  impactor  and  specimen.  In  the  next  step,  one  wave  reaches  the 
back  surface  of  the  impactor  while  the  other  wave  continues  into  the 
specimen.  The  wave  in  the  impactor  is  reflected  as  a rarefaction  wave 
and  propagates  back  through  the  impactor,  toward  the  specimen.  Sub- 
sequently, a rarefaction  wave  moves  into  the  specimen  from  its  free  sur- 
face. When  the  two  plane  rarefaction  waves  interact,  the  result  is  a 
tensile  stress  in  the  specimen.  This  tensile  stress  occurs  because  in 
rarefaction  waves,  material  is  accelerated  in  a direction  opposite  to 
the  rarefaction  propagation  direction.  If  the  tensile  stress  amplitude 
and  duration  are  sufficient,  incipient  spall  fracture  occurs. 


BEFORE  IMPACT 


Figure  7.  Simplified  schematic  of  the  spall  process.  After  impact 
the  time  sequence  of  wave  interactions  is  shown  from 
left  to  riqht.  The  solid  and  dashed  lines  indicate  the 
propagation  of  compression  and  rarefaction  waves, 
respectively . 

The  fracture  process  is  characterized  by  a threshold  stress 
for  crack  or  void  nucleation  as  well  as  the  other  material-dependent  NAG 
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parameters.  Thus,  controlled  stress  conditions  will  result  in  controlled 
amounts  of  spall-fracture  damage.  The  stress  amplitude  is  determined 
from  the  impact  velocity  and  the  previously  measured  shock  properties  of 
the  materials  used.  The  stress  duration  is  controlled  by  the  thickness 
of  the  impactor;  the  larger  the  thickness,  the  larger  the  duration  of 
the  stress. 

A summary  of  the  spall  fracture  experiments  performed  to  obtain 
the  NAG  parameters  for  HF-1  steel  is  presented  in  Table  3.  The  last 
column  of  the  table  gives  a qualitative  degree  of  damage  for  the  speci- 
mens as  determined  by  microscopic  observation.  Typical  tensile  stress 
amplitudes  ranged  from  about  2 to  5 GPa,  and  typical  stress  durations 
were  fractions  of  a microsecond. 


Table  3.  Spall  fracture  experiments  for  HF-1  steel. 


Number 

Projectile 

Velocity 

(km/s) 

89 

0.120 

95 

0.130 

111 

0.193 

112 

0.204 

110 

0.149 

92 

0.180 

113 

0.  202 

98 

0. 160 

91 

0.183 

90 

0.191 

94 

0.200 

88 

0.276 

Driver 

'hickness 

(mm) 

Sample 

Thickness 

(mm) 

1.609 

3.179 

2.369 

6.318 

1.156 

3.185 

1.159 

3.192 

2.371 

6.348 

2.372 

6.  336 

1.596 

3.189 

2.  372 

6.349 

2.373 

6.337 

2.370 

6.361 

2.369 

6.355 

1.596 

3.190 

Qualitative 
Degree  of 

Damage 

No  Visible  Damage 

Incipent 

Light 

Light 

Light 

Light 

Medium 

Medium 

Medium 

Heavy 

Heavy 

Full  Separation 


C.  PROCEDURE  FOR  OBTAINING  CRACK  DISTRIBUTIONS 


To  obtain  crack  distributions,  the  recovered  and  sectioned 
specimens  (Figure  8)  are  first  macrophotographed  at  5X  magnification 
using  a Bausch  and  Lomb  L camera.  During  this  procedure,  it  is  very 
important  to  illuminate  the  specimen  with  diffuse  light  so  that  all  the 
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cracks  are  visible  in  the  photograph.  Figures  9 and  10  are  composites  of 
these  photographs  for  selected  specimens.  For  crack  counting,  each  5X 
photograph  was  enlarged  again  and  arranged  together  to  form  a 25X  com- 
posite of  the  specimen  surface. 

SOFT  RECOVERED  SECTIONED 

Figure  8.  Sectioning  of  specimen  to  reveal  fracture  damage. 

The  25X  composite  was  mounted  on  a Computer  Equipment  Corpora- 
tion Model  TF-10C  Digitizer  to  determine  the  coordinates  of  the  cracks. 
The  general  data  reduction  criteria  used  was  to  count  the  longest  segment 
as  a crack  and  each  intersecting  segment  as  a crack.  If  the  intersecting 
crack  extended  through  the  primary  crack,  its  length  was  taken  as  extend- 
ing from  tip  to  tip.  If  the  intersecting  crack  stopped  at  the  primary 
crack  or  appeared  to  extend  through  the  primary  crack  with  a relative 
displacement,  each  segment  was  treated  as  an  individual  crack. 

The  digitized  data  were  then  used  as  input  to  a computer  pro- 
gram that  divides  the  damaged  region  into  ten  zones  parallel  to  the 
specimen  faces,  calculates  the  length,  orientation,  and  location  of  the 
midpoint  of  each  crack  relative  to  the  impact  surface,  and  assigns  each 
crack  to  a zone  based  on  the  location  of  the  midpoint.  This  program 
was  written  to  interface  with  the  SRI  BABS2  program.  The  computer- 
generated representation  of  the  digitized  surface  crack  data  for  a 
specimen  is  shown  in  Figure  11.  The  volume  distribution  of  cracks 
for  a specimen  is  shown  in  Figure  12.  This  distribution  is  represented 
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(a)  shot  ill,  0.193  km/s 


(b)  shot  112,  0.204  km/s 


(c)  SHOT  113,  0.202  km/s 


Figure  9.  Macrophotoqraphs  of  '.  2-tran- thick  specimens  showino 
microscopic  spall  fracture. 
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by  a plot  showing  the  number  of  cracks  having  a radius  R (crack  half- 
length)  greater  than  a certain  value.  To  eliminate  the  effects  due  to 
relief  waves  from  the  edge  of  the  specimens,  only  fracture  data  central 
to  the  damaged  region  were  used. 


I0*t 


Figure  11.  Computer  plot  of  the  surface  distribution  of  cracks  for 
a specimen.  The  marks  on  the  specimen  edge  indicate  the 
boundaries  for  the  ten  damage  zones. 

D.  CALCULATION  OF  NAG  PARAMETERS 

The  volume  crack  distributions  are  quantified  by  fitting  a 
straight  line  to  the  data  corresponding  to  zone  of  maximum  damage  for 
each  specimen.  This  straight  line  fit  (Figure  12)  is  represented  by  the 
equation 

N = N exp  (-R/R, ) (2) 

go  1 

where  is  the  number  of  cracks  per  unit  volume  with  radius  greater 

than  R,  Nq  is  the  total  number  of  cracks  generated  per  unit  volume,  and 

R,  is  related  to  the  final  size  of  the  cracks.  Values  of  N and  R,  are 
1 o 1 

determined  for  each  specimen. 

The  maximum  tensile  stress  and  tensile  duration  are  then  deter- 
mined for  each  specimen  using  the  PUFF  code  and  allowing  no  fracture  to 
occur.  The  quantities  obtained  from  the  computation  are  the  peak  tensile 

stress  o and  the  tensile  duration  At.  The  tensile  duration  is  the  time 
m 

for  which  a threshold  stress  value  is  exceeded. 
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R ( CRACK  RADIUS.  MM  1 


Figure  12.  Computer  plot  of  the  volume  distribution  of  cracks  for 
a specimen. 
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The  quantities  N^,  R^ , At,  and  are  used  to  form  the  first  estimates 
of  the  NAG  parameters.  The  nucleation  rate  for  cracks  N and  the  growth 
rate  for  cracks  R are  given  by 


fa  - a ~i 

k = No  e4  ”ain0J 


R = T,  (0-0  ) R (4) 

1 go 

where  a is  the  applied  stress  and  R is  the  crack  radius.  The  NAG  param- 
eters are  N , a, . a , o , T, , and  R . R is  the  crack  radius  at  nuclea- 
o 1 no  go  1 o o 

tion  and  is  defined  in  Equation  (5). 

A first  estimate  for  the  threshold  nucleation  rate  N and  the 

o 

nucleation  sensitivity  o^  is  obtained  by  plotting  the  values  lnlN^/At) 

versus  o^  for  each  of  the  specimens  and  fitting  the  points  to  Equation 

(3) . o is  the  threshold  stress.  The  growth  coefficient  T and  the 
no  1 

nucleation  size  parameter  Rq  are  obtained  by  plotting  the  values  InRj, 

versus  o At  for  each  of  the  specimens  and  fitting  the  points  to  the 
m 

equation 


InR,  = InR  + T,  o At 
1 o 1 m 


which  is  obtained  by  integrating  Equation  (4) , assuming  the  growth 

threshold  parameter  to  be  negligible.  The  material  viscosity  i|  is 

obtained  from  the  equation  T = l/4g.  The  dynamic  plain-strain  fracture 

1 1/2 

toughness  Kt  , is  obtained  from  the  equation  n = ( tt/4R)  K_. 

IC  qo  IC 

These  NAG  parameter  estimates  are  then  used  as  input  to  a PUFF 

computation,  which  allows  fracture  to  occur.  New  calculated  values  for 

N and  R,  are  obtained  from  the  PUFF  output  and  compared  with  the  expert  - 
o 1 

mental  values.  If  agreement  is  not  obtained,  the  parameters  N and  T 


are  adjusted  and  the  computation  is  repeated.  This  iterative  process 
is  continued  until  the  experimental  and  computed  values  agree.  The  NAG 
parameters  determined  for  HF-1  steel  are  given  in  Table  4.  Also  included 
for  comparison  are  the  SRI-determined  parameters  for  Armco  iron. 


Table  4.  Dynamic  fracture  parameters  for  HF-1  steel  and  Armco  iron. 


Parameter 


Units 

cm2/dyn-s 

dyn/cm2 

cm 

no. /cm^-s 

dyn/cm2 

dyn/cm2 


HF-1  Steel 


5.25  x 10 
1.0  x 109 


5.5  x 10 

4.0  x 10 

2.0  x 10 


-3 

10 


10 


1.0  x 10' 


Armco  Iron* 

6.0  x 10"4 

2.0  x 108 

5.0  x 10-5 
4.6  x 1012 

3.0  x 109 
4.56  x 109 


* Reference  1 

V.  EXPLODING  CYLINDER  EXPERIMENTS 
A.  INTRODUCTION 


Two  types  of  experiments  were  performed  to  provide  information  on 
the  expansion  and  fragmentation  of  explosive-filled  cylinders:  (1)  saw- 

dust pit  experiments  to  provide  fragments  exhibiting  both  radial  and  shear 
fractures  for  model  development  and  to  provide  the  fragment  mass  distri- 
butions, and  (2)  a framing  camera  experiment  to  provide  cylinder  strain 
data.  The  test  configurations  consisted  of  cast-in-place  composition  B 
explosive  which  extended  beyond  the  cylinder  to  minimize  end  effects 
(Figure  13) . Configuration  details  are  given  in  Table  5. 
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Figure  13. 


Schematic  of  explod 


Table  5.  Configuration  details 


Experiment  Test 


Number 

Type 

Material 

(cm) 

(cm) 

Charge 

1 

Pit 

Armco  Iron 

7.62 

11.43 

20.  3 

2.77 

2 

Pit 

Armco  Iron 

7.62 

11.43 

20.  3 

2.76 

3 

Pit 

HF-1  Steel 

7.62 

12.07 

20.  3 

2.75 

4 

Pit 

HF-1  Steel 

7.62 

12.07 

20.  3 

2.75 

5 

Framing 

Armco  Iron 

7.62 

11.43 

38.1 

4.10 

Camera 
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B.  SAWDUST  PIT  EXPERIMENTS 

These  experiments  were  conducted  at  the  NSWC/DL  Sawdust  Pit 
Facility.'*'*’  Following  detonation,  the  fragments  were  captured  in  a 
large  volume  of  sawdust,  and  then  were  recovered  by  screen  and  magnetic 
separation  techniques.  Mass  distributions  were  determined  by  weighing 
and  counting  each  fragment  with  mass  greater  than  1 grain  (64.8  mg) . 
Fragments  smaller  than  this  were  not  counted,  but  the  total  weight  was 
recorded.  In  all  experiments,  at  least  99.5  percent  metal  mass  was 
recovered. 

Fragment  mass  distributions  determined  from  the  pit  experiments 
are  shown  in  Figures  14  and  15.  Both  the  Armco  iron  and  HF-1  steel 
distributions  deviated  from  linear  behavior  at  large  fragment  masses. 

This  is  not  unusual  because  of  the  small  statistical  sample  in  the 
highest  mass  classes.  The  Armco  iron  cylinders  also  deviated  from 
linear  behavior  in  the  low  mass  classes.  This  is  considered  to  be  an 
actual  feature  of  the  distribution  and  not  due  to  secondary  breakup  of 
the  fragments  during  recovery  since  Armco  iron  is  a tough  and  ductile 
material  compared  with  the  HF-1  steel  (which  did  not  deviate  from 
linear  behavior) . 

C.  METALLOGRAPHIC  OBSERVATIONS  OF  RECOVERED  FRAGMENTS 

The  microstructural  features  of  representative  fragments  were 
investigated  by  metallographic  procedures.  Selected  fragments  of  Armco 
iron  and  HF-1  steel  are  shown  in  Figure  16.  The  fragments  had  well- 
defined  fracture  surfaces  which  were  preserved  with  Krylon  lacquer. 
Laboratory  specimens  were  prepared  by  sectioning  the  fragments  to 
reveal  a through-thickness  plane  normal  to  the  cylinder  axis.  The 
specimens  were  mounted  and  polished,  and  then  etched  (2-percent  Nital 
solution)  to  reveal  the  microstructural  features  of  interest.  Specimen 
examination  was  performed  with  a Bausch  and  Lomb  Research  II  Metal lograph. 
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NSWC/DL  (GFfj) 

i — 3 INCHES -■* 


Figure  10>.  Fragments  of  Armco  iron  (upper  in  fiqure)  arid  IIF-1 
steel  (lower  in  fiqure)  used  in  metal lographic  and 
SUM  studies. 

Fracture  surfaces  were  examined  with  a Cambridge  MARK  I I -A 
scanning  electron  microscope  (SEM).  To  prepare  the  SI'M  specimens,  the 
clear  lacquer  was  removed  by  soakinq  the  fragments  in  acetone.  The 
fracture  surfaces  were  then  cleaned  ultrasonie.il  ly  in  a solution  consist- 
ing of  1-ml  hydrochloric  acid  (sp.  qr.  1.19),  4 -ml  .’-Butyno- 1 , 4-diol 
(1.1  percent  aqueous  solution),  and  10-ml  deionized  water.  This  method 
also  removed  cellulose  fibers  from  the  sawdust  which  had  become  imbedded 
in  the  fracture  surfaces  during  recovery.  The  fragments  wore  cemented  to 
aluminum  disks  with  an  electrically-conducti ve  silver  paste  so  that  the 
fracture  surfaces  could  bo  viewed  in  the  SKM. 

Fiqure  17  shows  cross-section  specimens  of  Armco  iron  ami  III-'—  L 
steel.  Both  specimens  exhibited  bi modal  fracture.  For  convenience  of 
discussion,  the  specimens  are  divided  into  three  regions  in  a radial 

In  Region  1 (approximately  the  outer  one-half  thickness)  the 


direction. 


specimens  exhibited  transgranular  cleavage  on  planes  roughly  normal  to  the 
tangential  direction.  While  both  metals  showed  internal  cracking,  the 
cracks  in  HF-1  steel  were  fewer  in  number  and  longer  in  length.  These 
cracks  remained  sharp,  with  no  evidence  of  crack-tip  blunting  due  to 
plastic  deformation  [Figure  18(b)].  In  Armco  iron  the  cleavage  cracks 
tended  to  be  stacked  in  radial  arrays  with  each  individual  crack  approxi- 
mately one  grain  diameter  in  length  [Figure  18(a)].  The  cracks  exhibited 
blunting  which  progressed  as  the  location  of  the  cracks  approached  the 
transition  Region  2.  Those  cracks  near  Region  2 were  highly  lenticular, 
suggesting  a strong  plastic  deformation  gradient. 


(a) 


(b) 


Figure  17.  Cross-sections  of  recovered  fragments:  (a)  Armco  iron 

and  (b)  HF-1  steel. 


i 
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Figure  18.  Internal  cracks  in  Region  1:  (a)  Armco  iron  and 

and  (b)  HF-1  steel. 

Figure  19(a)  is  an  SEM  micrograph  of  a radial  fracture  surface  of 
Armco  iron  in  Region  1.  The  fracture  morophology  consists  of  regions 
of  quasi-cleavage  and  pare  shear.  This  observation,  in  conjunction  with 
Figure  18(a),  suggests  that  in  Armco  iron  the  crack  nucleation  and  the 
crack  growth  processes  are  different.  Nucleation  apparently  occurs  by 
cleavage  of  individual  grains,  whereas  crack  growth  and  link-up  occur  by 
void  nucleation,  growth,  and  coalescence.  Figure  19(b)  is  an  SEM  micro- 
graph of  a radial  fracture  surface  in  HF-1  steel.  The  fracture  morphol- 
ogy again  consists  of  regions  of  quasi-cleavage  and  regions  of  void  coa- 
lescence. This  morphology  compares  favorably  with  that  observed  in 
spalled  gas-gun  specimens.  The  scale  of  the  fracture  is  finer  in  the 
fragments  than  in  the  spall  specimens,  and  the  spall  specimens  exhibit 
more  cleavage.  The  difference  in  the  two  cases  may  result  from  the 
larger  hydrostatic  compression  that  occurs  during  cylinder  fracture. 
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Seaman  et  al.  have  observed  crack  blunting  in  Armco  iron  frag- 
ments supplied  to  them  for  metallurgical  examination.  They  suggest  that 
this  fracture  morphology  is  produced  by  tensile  cleavage  of  individual 
grains,  the  arrest  of  these  cracks  at  grain  boundaries,  and  the  subse- 
quent plastic  stretching  at  the  crack  tips. 


Figure  19.  SEM  micrographs  of  radial  fractures  in  Region  1: 
(a)  Armco  iron  and  (b)  HF-1  steel. 


Another  important  feature  of  the  microstructure  in  Region  1 is 
the  distortion  of  the  matrix  material  between  radial  fractures.  Micro- 
hardness measurements  across  the  thickness  (Figure  20)  indicate  the 
relative  degrees  of  material  hardening.  Armco  iron  shows  a strong  strain 
gradient  in  the  middle  third  of  the  thickness,  whereas  HF-1  steel  shows 
a more  gradual  trend.  To  estimate  the  tangential  strain  in  this  reqion, 
one  can  consider  constant-volume  spherical  grains  shearing  to  form 
ellipsoids  of  the  same  volume.  If  x and  y are  the  major  and  minor 
ellipsoid  axes,  respectively,  the  tangential  strain  is  approximately 
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Figure  20.  Microhardness  profiles  through  the  fragment  thickness. 

Since  the  virgin  Armco  iron  (Figure  3)  had  an  equiaxed-grain 

structure,  Equation  (6)  can  be  used.  An  average  over  15  grain  diameters 

gives  e =0.33.  Since  stress  relief  associated  with  the  formation  of 
x 

these  fractures  arrests  matrix  deformation,  it  is  estimated  that  the 

ruuial  fractures  were  formed  approximately  13  gs  after  detonation. 

19 

Bardes  obtained  f = 1.8  for  an  HF-1  steel  cylinder  with  a carbide 
network  and  well-tempered  martensite  matrix. 

In  Region  3,  the  fracture  surfaces  were  of  the  shear  type  in 
both  materials.  The  Armco  iron  specimens  (Figure  21)  exhibited  larqe 
grain  distortions  in  the  bulk,  and  unstable  shear  at  the  fracture  surfaces 
and  also  at  intervals  within  the  bulk.  Matrix  deformation  occurred  by 
twinning  and  glide  on  closely  spaced,  wavy  slip  bands  within  each  grain. 
Regions  of  highly  localized  shear  were  inclined  approximately  45°  to  the 
cylinder  radius,  following  the  maximum  shear  trajectories  in  discrete 
packets.  The  deformation  became  less  localized  at  the  ends  of  these 
packets.  At  the  intersection  of  each  active  shear  trajectory  with  the 
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inner  surface,  a finite  offset  was  observed.  Seaman  et  al.  ‘ , however, 
reported  instances  of  active  shear  trajectories  with  no  apparent  offset 
at  the  inner  surface. 

RADIAL 

j DIRECTION 


Figure  21.  Micrograph  showinc  localized  shear  deformation  in 
Armco  iron. 

The  spacing  between  the  active  shear  trajectories  in  Armco  iron 
that  are  bounded  by  fracture  surfaces  of  the  same  rotation  was  approxi- 
mately 1 to  5 mm.  The  shear  packet  width  was  approximately  200  pm,  and 
the  offset  at  the  inner  surface  was  typically  390  urn. 

Observations  on  HF-1  fragments  in  Region  3 were  somewhat  different: 
bulk  distortions  were  much  less  severe  and  localized  shoarinq  was  of  the 
transformation  type.  This  phenomena,  which  lias  been  widely  observed  in 
the  high  strain  rate  processing  of  steels,  has  been  explained  as  localized 
shear  heating  to  temperatures  in  the  austenite  phase  field.  When  the 
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deformation  is  complete  (10  to  50  ;is  after  detonation) , there  is  a rapid 
quenching  of  the  heated  zone  by  the  surrounding  cooler  material.  This 
forms  untempered  martensite  which  is  not  attacked  by  the  etching  material 
because  of  the  lack  of  c-carbide  precipitates  in  the  substructure.  Hence, 
a large  portion  of  the  shear  zone  is  clearly  delineated.  Figure  22(a) 
shows  a portion  of  a shear  zone  in  HF-1  steel,  and  a crack  associated 
with  the  untempered  martensite. 

Along  a single  shear  trajectory  which  has  a crack  associated  with 
it,  there  are  often  regions  exhibiting  crack  openings  in  the  bulk.  This 
observation  for  HF-1  steel  [Figure  22(b)]  corresponds  to  the  stop-start 
behavior  in  Armco  iron  and  again  suggests  interior  nucleation. 

Shear  zones  of  different  rotation  sometime  intersect  as  shown  in 
Figure  22(c).  When  this  occurs,  a convenient  marker  is  provided  for 
offset  measurements.  Here  the  offset  associated  with  the  primary  shear 
is  250  im,  corresponding  to  an  offset  of  350  pm  at  the  inner  surface 
of  the  same  primary  shear.  The  spacing  between  active  shear  zones  of 
the  same  rotation  was  between  320  and  1450  urn  and  the  shear  packet  width 
in  uncracked  shear  zones  was  2.5  to  6.8  pm.  The  above  metal lographic 
observations  are  important  for  computational  models  of  the  shear  phenomena 
and  will  be  discussed  further  in  the  next  chapter. 

An  SEM  microqraph  of  the  shear  fracture  in  Armco  iron  is  shown  in 
Figure  23(a),  where  the  fracture  has  occurred  by  void  growth  and 
coalescence.  Figure  23(b)  shows  a higher  magnification  of  one  of  the 
spheres  in  Figure  23(a).  Quantitative  dispersive  X-ray  analysis  of  the 
spheres  failed  to  reveal  any  difference  in'  compos i t ion  between  this 
feature  and  the  matrix  material.  It  is  concluded  that  t he...  feat  ares 
resulted  from  the  col  lapse  of  very  hot  asperr  i follewin:  • s it  it  ion. 

An  SEM  microqraph  of  the  shear  fracture  in  HF-1  stool  i-  h>  -wti  in  ! i .!■:•■  .'-j . 
The  fracture  again  resulted  from  void  coalescence.  In  'hi  , >weve»  , 

f he  scale  is  much  finer  Mian  in  Armco  iron. 
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Figure  23.  SEM  micrographs  of  the  shear  fracture  surface  from  an 
Armco  iron  fragment  showing:  (a)  fracture  by  void 

growth  and  coalescence  and  (b)  high-magnification  view 
of  melted  asoerrites. 
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Figure  24.  SUM  micrograph  of  the  shear  fracture  surface  of  an  UK-1 
steel  fragment. 
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D.  FRAMING  CAMERA  EXPERIMENT  AND  HEMP  CODE  PREDICTIONS 


An  attempt  was  made  to  measure  the  strain  history  of  an  expand- 
ing cylinder  using  the  Moire  interference  technique.  The  38.1-cm-long 
Armco  iron  cylinder  (Table  5)  was  used  in  this  test.  Prior  to  filling 
the  cylinder  with  explosive,  a grid  pattern  and  a parallel  line  pattern 
were  silk  screened  onto  the  surface  as  shown  in  Figure  25.  The  large 
grid  pattern  consisted  of  1.6-mm-wide  lines  over  a square  coordinate 
grid  of  6.35  mm.  The  grid  extended  6.35  cm  along  the  cylinder  axis  and 
17.8  cm  around  the  circumference.  The  parallel  line  pattern  consisted  of 
a series  of  0.5-mm-wide  lines  12.7-cm  long,  oriented  along  the  axis  of 
the  cylinder.  The  line  spacing  was  1 mm  and  the  pattern  also  extended 
17.8  cm  around  the  circumference.  The  patterns  were  centered  along  the 
axis  of  the  cylinder. 


Figure  25.  HF-1  steel  cylinder  with  grids  for  strain  measurements. 

A schematic  for  investigating  interference  patterns  is  shown  in 
Figure  26.  A photograph  of  the  initial  pattern  was  used  as  the 
reference  grating.  This  arrangement,  devised  by  SRI,  provides  Moire 
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fringes  which  can  be  directly  recorded  by  the  framing  camera  when  focused 
on  the  plane  of  the  analyzing  grid.  It  was  determined,  however,  that  this 
procedure  could  not  be  used  because  the  focusing  of  the  framing  camera  at 
the  analyzer  plane  required  the  camera  to  be  too  close  to  the  exploding 
cylinder. 


HIGH  RESOLUTION  ANALYZER  AND  PLANE 

LENS  OF  INTERFERENCE  PATTERN 


Figure  26.  Schematic  for  optical  measurement  of  strains 
(Reference  4) . 


The  procedure  that  was  used  for  the  strain  test  was  as 
follows:  Before  detonation,  a series  of  analyzer  patterns  was  made  with 

the  framing  camera.  On  detonation,  the  event  was  recorded  at  a 4.17-iis 
per  frame  rate.  The  analyzer  patterns  were  then  super imposed  on  the 
detonation  photographs.  This  procedure  should  have  produced  Moire  inter- 
ference fringes,  but  there  was  insufficient  detail  in  the  film  to  resolve 
the  patterns.  It  was,  nevertheless,  possible  to  obtain  the  gross  strain 
history  of  the  cylinder  by  measuring  its  diameter  on  the  framing  camera 
photographs,  such  as  Figure  27. 

HEMP  code  calculations  were  performed  for  the  expansion  of  both 

Armco  iron  and  HF-1  steel  cylinders.  The  results  are  given  in  Figure;? 

2f> 

2H  and  20,  respectively.  Previous  work  by  Coughlin  has  shown  that 

HEMP  code  calculations  predict  the  expansion  of  a cylinder  under  exj  losive 

loading  to  within  about  10  percent.  To  minimize  comput.it ional  difficulties, 
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Figure  27.  Framing  camera  photograph  of  Armco  iron  cylinder,  56  ;.s 
after  detonation. 

the  point  of  initiation  was  placed  at  the  end  of  the  metal  cylinder. 

Strain  effects  on  the  yield  stress  for  Armco  iron  were  approximated  by 

assuming  the  yield  stress  increased  from  180  to  350  MPa.  This  assumption 
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is  quite  reasonable  as  shown  by  the  data  of  Hocket  and  Zukas.  Figure  30 
is  a comparison  of  the  HEMP  code  calculations  with  results  of  the  framing 
camera  experiment  for  Armco  iron. 

VI.  SHEAR  BAND  MODELS  AND  FRAGMENT  MASS  DISTRIBUTION  COMPARISONS 
A.  INTRODUCTION 

In  this  chapter,  two  models  that  have  been  developed  for 
shear  band  formation  in  cylinder  fragment  at  ion  are  discussed.  They  art- 
based  in  part  on  t lie  meta 1 lographie  observat  ions  that  were  presented  in 
the  pre  ceding  chapter.  fn  addition,  eonipar i sons  are  presented  for 
measured  and  computed  fragment  mass  distributions  foi  explosive-filled 
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OUTSIDE  RADIUS  AT  ONE-HALF  CYLINDER  HEIGHT,  CM 

Figure  30.  Comparison  of  framing  camera  results  and  HEMP  code 
predictions  for  Armco  iron  cylinders. 

cylinders  of  Armco  iron  and  HF-1  steel.  The  computed  mass  distributions 

were  obtained  by  SRI  using  their  recently  developed  model  for  cylinder 
4 

fragmentation. 

B.  SHEAR  BAND  MODELS 

Highly  localized  shear  deformation  has  frequently  been  observed 

in  metals  deformed  at  high  strain  rates'^  or  at  cryogenic  temperatures.^'^3 

This  type  of  deformation  is  often  referred  to  as  adiabatic  shear.  Zener 
31 

and  Holloman  explain  the  phenomena  as  an  instability  induct'd  bv  the  heat 
generated  during  plastic  flow.  The  heat  produces  sufficient  thermal  soft- 
ening to  offset  the  effects  of  strain  and  strain-rate  hurdenina  anti  onuses 

the  deformation  to  concentrate  in  narrow  bands.  There  are  several  analytical 

3',-go 

treatments  of  the  problem  in  the  literature;  however,  none  of  them 

contain  the  required  kinetic  information  necessary  for  calculatin'!  the  pro- 
longation of  shear  instabilities. 
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A dislocation  theory  of  unstable  shear  deformation  has  been 
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formulated  by  Arsenault  and  Crowe.  This  theory  is  based  on  the  obser- 
vation that  plastic  deformation  in  many  metals  and  alloys  is  thermally 
activated.  Consider  a small  cylindrical  mass  of  material  with  a disloca- 
tion source  on  a slip  plane  P as  shown  schematically  in  Figure  31,  the 
shear  strain  rate  ^ is  given  by  the  Orowan  equation 

Y = apbv  (7) 

where  a is  the  geometric  factor  taken  as  1/2,  b is  the  magnitude  of  the 
Burgers  vector,  p is  the  density  of  mobile  dislocations  on  the  slip 
plane,  and  v is  the  average  dislocation  velocity.  p can  be  written  as 
the  mobile-line  length  per  unit  volume  according  to 


N ( y ) tt  r 
V 

o 


(8) 


where  N(y)  is  the  number  of  dislocations  (a  function  of  strain),  r is  the 

dislocation  loop  radius,  and  V is  the  unit  volume.  If  the  deformation 

o 

proceeds  by  nucleation  and  expansion  of  dislocation  loops  from  small 
nucleate  radii  to  larger  values  r,  then  r = vt,  where  t is  the  time  of 
the  loop  expansion.  Then, 


Figure  31. 


Schematic  of  shear  deformation  via  dislocation  motion 
on  a slip  plane. 


40 


According  to  the  theory  of  thermally-activated  plastic  deforma- 
42 

tion,  the  dislocation  velocity  is  given  by  an  Arrhenius  equation  of 
the  form 


(-AH<T,o)\ 
XT  ) 


where  AH  is  the  activation  enthalpy  which  is  a function  of  the  stress 

s and  the  absolute  temperature  T,  k is  the  Boltzmann  constant,  and  v is 

o 

the  maximum  dislocation  velocity. 

Considering  dislocation  motion  as  a power  dissipation  process, 
the  power  dissipated  in  a small  annular  region  about  the  dislocation  line 


P = oyV 


where  V is  the  volume  of  the  annular  region  given  approximately  by 

2nr6b,  where  <Sb  is  the  cross-sectional  area  of  the  annulus.  The  temperature 

43 

distribution  resulting  from  the  radiating  line  source  is  given  as  a 
function  of  position  (x,y)  by: 


T (x , y) 


7R 


where  K^lr/X)  is  a modified  Bessel  function,  \ = 2H/V  is  a characteristic 
distance,  H = k/C  is  the  thermal  diffusivity,  is  the  thermal  conduc- 
tivity, and  C is  the  specific  heat  at  constant  volume.  Since  the  dis- 
v 

location  velocity  increases  with  increasing  temperature,  the  passaqe  of 
the  first  dislocation  heats  the  material  on  the  slip  plane.  The  second 
dislocation  emitted  from  the  source  thus  travels  at  a slightly  hiqher 
velocity  than  the  first  dislocation.  If  the  strain  rate  is  sufficiently 
high,  this  effect  can  cascade  into  unstable  shear  deformation  upon  sub- 
sequent dislocation  emission  from  the  source. 
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The  above  equations  were  programmed  for  an  iterative  solution 

using  a finite-difference  technique  to  solve  the  heat  flow  equations. 

The  calculations  were  not  completed,  however,  because  the  computer  time 

required  for  a solution  was  prohibitively  long.  In  principle,  the  above 

concepts  provide  a physically  realistic  description  of  the  onset  and 

propagation  of  unstable  thermoplastic  shear  deformation.  Campbell 
44 

et  al.  have  shown  that  dislocation  sources  begin  to  operate  in  about 
5 ns.  Source  operation  is,  therefore,  probably  not  the  rate-controlling 
factor,  provided  the  dislocations  leaving  the  source  do  not  apply  suffi- 
cient back  stress  to  turn  off  source  operation.  For  pressure  durations 
associated  with  explosive  loading,  this  is  not  expected  to  occur. 

The  above  model  is  also  consistent  with  the  metallographic 
observations  described  earlier  in  this  report.  In  particular,  the  model 
provides  a mechanism  for  the  nucleation  of  shear  zones  within  the  bulk 
material,  since  shear  nuclei  will  occur  in  grains  having  slip  systems 
favorably  oriented  with  the  maximum  shear  trajectories.  The  difference 
between  the  SRI  and  NSWC/DL  observations  regarding  offset  at  the  inner 
surface  are  not  significant  since  they  may  only  reflect  the  proximity  of 
favorably  oriented  grains  to  this  surface. 

An  important  question  that  arises  concerns  the  initiation  of 
an  instability  at  any  particular  location.  One  explanation  is  that  each 
zone  is  associated  with  a local  inhomogeniety  which  intensifies  the  stress 
on  that  trajectory.  Another  explanation  is  that  those  trajectories 
containing  the  largest  number  of  grains  favorably  oriented  for  slip  are 
naturally  weaker,  and  become  preferred  sites  of  unstable  shear. 
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The  SRI  approach  ' to  the  localized  shearing  problem  com- 
bines the  bulk  thermodynamic  approach  and  the  dislocation  approach.  This 
model  uses  a geometrically-scaled  analogue  of  the  dislocation.  Figure  12 
is  a schematic' of  the  model.  W is  the  width  of  the  shear  band  and  B is 
the  offset  associated  with  shear  band  propagation.  The  kinetics  of  pro- 
pagation are  again  described  by  the  Orowan  equation 
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where  v is  the  shear  strain  rate,  N , is  the  density  of  mobile  shear 

sb 

bands  and  v , is  the  averaqe  shear  band  velocity.  In  this  formalism  the 

. Sb  . 45 

governing  equations  are 


(14) 


(15) 


(16) 


where  p is  the  density,  Cv  is  the  specific  heat  at  constant  volume,  k 
is  the  thermal  conductivity,  Y is  the  yield  stress,  AE^  = is  the 

plastic  work,  and  E is  the  internal  energy  required  to  remove  most  of 
the  material  strength. 
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Figure  1.’ . 


Schema  tic  of  shear  instability  us  in  cal 
dislocation  motion  (Reference  45). 


Values  of  W and  B determined  metallographically  from  recovered 

fragments  are  given  in  the  previous  chapter.  These  quantities  should 

be  obtainable  by  performing  small-scale  shearing  experiments  such  as 
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those  described  by  Culver  [i.e.,  calculable  from  mechanical  property 
and  thermodynamic  data  and  Equations  <14),  (15)  and  (16)].  At  present, 

very  little  work  of  this  type  has  been  performed. 

Another  important  feature  of  these  models  is  that  the  shearing 
phenomena  is  primarily  associated  with  a plastic  deformation  process  and 
not  a fracture  process.  In  principle,  the  shear  mechanism  can  operate 
to  form  new  surfaces  in  perfect  Mode  II  fracture  propagation  (the  shear 
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component  of  stress  is  applied  normal  to  the  leading  edge  of  the  crack) . 
The  SEM  micrographs  of  Figure  23  and  24  indicate  that  this  phenomena 
does  not  occur  in  Arm co  iron  and  HF-1  steel.  In  these  materials,  the 
separation  associated  with  the  shearing  phenomena  probably  occurs  during 
the  later  stages  of  the  fragmentation  process  as  the  tangential  stress 
becomes  tensile  in  the  inner  half-thick.. ess  regions. 

C.  COMPARISON  OF  COMPUTED  AND  MEASURED  FRAGMENT  MASS  DISTRIBUTIONS 

In  this  section,  the  SRI  f ragmentation  model  computations 
are  compared  with  the  NSWC/DL  sawdust  pit  experiments.  The  procedure 
used  to  compute  a fragment  mass  distribution  for  an  exploding  cylinder  is 
as  follows:** 

(1)  Use  PUFF  to  calculate  the  elastic-plastic  stress  state  in 
each  undamaged  element  of  mass. 

(2)  At  each  time  step,  compare  computations  with  criteria  for 
brittle  fracture  and  shear  band  formation. 

(3)  If  the  shear  band  criteria  are  met  first  in  an  undam, .god 
mass  element,  use  the  shear  band  model  for  subsequent  computations. 

(4)  If  the  brittle  fracture  criteria  are  met.  first  in  an 
undamaged  mass  element,  use  the  brittle  fracture  model  for  subsequent 
eomput a f ions . 

(5)  Construct  the  fragment  mass  distribution  by  add i n i the 
crack  distributions  for  all  cells  and  weighting  the  sum  by  the  cell  mas.  . 
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The  comparison  of  the  mejsure'l  and  calculated  fragment  mass 
distributions  for  Armco  iron  and  HF-1  steel  are  presented  in  Figures  33 
and  34,  respectively.  As  can  be  seen  from  the  figures,  the  brittle 
fracture  and  shear  band  models  have  provided  reasonable  predictions  of 
the  fragment  mass  distributions  for  these  materials. 
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Comparison  of  measured  and  SR  I -ou  1 -ti  1 a t . fragment 
mass  distribution::  for  Arm  n iron  (Reference  4). 


Figure  13. 


VII.  SUMMARY  AND  RECOMMENDATIONS 


The  results  of  this  joint  NSWC/DL  - SRI  experimental  and  analyti- 
cal program  demonstrate  that  the  SRI  fragmentation  models  can  be  used  to 
calculate  fragment  mass  distributions  that  agree  favorably  with  experi- 
ment. These  results  verify  the  capability  of  newly  formulated  fracture 
and  shear  band  models  to  describe  fragmentation  behavior.  Material 
parameters  needed  in  the  models  are  obtained  from  dynamic  experiments 
and  thermodynamic  data. 

Further  model  refinement  is  needed  to  better  match  the  experimental 
results.  Additional  experimental  work  needed  includes  measurement  of  the 
NAG  parameters  for  other  heat  treatments  to  determine  the  sensitivity  of 
the  fragmentation  model  to  changes  in  material  microstructure  and  com- 
parison of  the  calculated  fragment  mass  distributions  for  these  heat 
treatments  with  sawdust  pit  data. 
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